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^ ! 1. Introduction 

o , 

, A Schrodinger operator is a differential operator whose symbol is the Laplace's 

operator. A quantum integral of a Schrodinger operator is a differential operator 
that commutes with it. 
p • A Schrodinger operator in m variables is called integrable if it has m algebraically 

D . independent quantum integrals in involution (i.e. commuting with each other). This 
^ ' notion is the quantum analogue of the notion of Liouville integrability of a classical 
Hamiltonian system. 

^ ■ A Schrodinger operator in m variables is called algebraically integrable if it 

■ is integrable but the algebra of its quantum integrals cannot be generated by m 
operators. In the one-variable case, algebraically integrable operators correspond 
to finite-gap potentials [Kr]. 

One of the most interesting examples of an integrable Schrodinger operator is the 
Calogero-Sutherland operator [C] , [S] , [OP] . This is the Hamiltonian of the quantum 
many-body problem with rational, trigonometric, or elliptic interaction potential. 
The Calogero-Sutherland operator depends on a parameter which is called the 
coupling constant. 

It has been observed [CV1],[CV2],[VSC] that the Calogero-Sutherland operators 
become algebraically integrable when the coupling constant takes a discrete set 
of special values. This is proved for the rational and trigonometric case but still 
remains a conjecture in the elliptic case for two or more variables. 

These results can be generalized to Calogero-Sutherland operators associated 
with root systems, which were defined in [OP]. 

In this paper we study integrability and algebraic integrability properties of cer- 
tain matrix Schrodinger operators. More specifically, we associate such an operator 
(with rational, trigonometric, or elliptic coefficients) to every simple Lie algebra 
g and every representation U of this algebra with a nonzero but finite dimen- 
sional zero weight subspace. (The Calogero-Sutherland operator is a special case 



of this construction). Such an operator is always integrable [E]. Our main result is 
that it is also algebraically integrable in the rational and trigonometric case if the 
representation U is highest weight. This generalizes the corresponding result for 
Calogero- Sutherland operators ([CVl]). We also conjecture that this is true for the 
elliptic case as well, which is a generalization of the corresponding conjecture from 
[CV2]. 

The proof of the main result is based on the method of i/^-function - a joint 
eigenfunction of quantum integrals of the Schrodinger operator. This method was 
developed in [CVl]. The proof of existense and uniqueness of the •(/^-function is 
based on an explicit construction of this function which uses representation theory 
of the Lie algebra g. To be more precise, the V'-function is realized (up to a factor) 
as a weighted trace of an intertwining operator between a Verma module over g and 
the tensor product of this module with U. Such realization goes back to [E],[EK1], 
where it is found that joint eigenfunctions of quantum integrals of a Calogero- 
Sutherland operator can be realized as traces. Using the theory of Shapovalov form 
for g ([Sh],[KK]), we prove that the trace function satisfies the axioms for the t/^- 
function analogous to those formulated in [CVl], and is determined uniquely by 
them, and then establish algebraic integrability using the method of [CVl]. 

The paper is organized as follows. In Section 2 we make the necessary definitions, 
motivate them, and formulate the main result. In Section 3 we give information 
about Verma modules, Shapovalov form, and intertwining operators. In Section 4 
we define the i/^-function as a normalized trace, and prove two properties of this 
function. In Section 5 we prove that these two properties uniquely determine the 
V'-function. In Section 6 we prove algebraic integrability using the V'-function. In 
the Appendix we describe how to get Weyl group invariant quantum integrals from 
central elements (Casimirs) of U (g) . 

Acknowledgements. The authors would like to thank I.Frenkel, A.Kirillov 
Jr., F.Malikov and A.Veselov for numerous discussions and useful suggestions. 

The first author's work was supported by the Alfred P.Sloan Graduate Disser- 
tation Fellowship. 

2. Main definitions and results 

Let y be a finite-dimensional complex vector space. 

Definition 2.1. A matrix differential operator is a differential operator whose coef- 
ficients are EndV -valued functions. A matrix Schrodinger operator is a differential 
operator of the form 

(2-1) L = A-yl(a;), a; e C^, 

where A is the Laplacian in C"^, and A is a meromorphic function in with 
values in EndV. 

Definition 2.2. A matrix Schrodinger operator L is called integrable if there exist 

pairwise commutative matrix differential operators Li = L, L2, such that the 

symbols of Li have the form pi{^^, -^—)Id, where pi are algebraically indepen- 
dent polynomials (pi{y) — y^). The operators Li,...,Ljn are called the quantum 
integrals for L. 

Let R G i)* — C"^ be the root system of a simple Lie algebra g of rank m, and 
let A+ be the set of positive roots of R. 
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Definition 2.3. The Calogero- Sutherland ( CS) operator for R is the operator 



(2-2) L = A- ^ Caw(< a;, a >) + 

aeA+ 

where the scalar constant may depend only on the length of the root a, and 
u is one of the following potential functions: (i) u{x) = 2jx^ (rational potential), 
(a) u{x) — 2/sinh'^x (trigonometric potential), or (Hi) u{x) — 2p{x\u>i,u>2) (elliptic 
potential), where p{x\ui, U2) is the Weierstrass elliptic function with periods ui, uj2, 
and K is a constant (Cases (i) and (ii) are degenerations of case (Hi)). 

Such operators were introduced by Calogero [C] and Sutherland [S] for the root 
system Am and by Olshanetsky and Perelomov [OP] in general. 

Theorem 2.1. The operator L given by (2-2) is integrable. Furthermore, the 
symbols of the quantum integrals Li, 1 < i < m, generate the algebra of Weyl group 
invariant polynomials on i). 

For R = Am (and in some other special cases) this theorem was proved in [OP] . 
Cases (i) and (ii) for general root systems were settled by Heckman and Opdam 
[HO, HI, 01, 02]. Case (iii) for Bm, Cm and Dm is settled in [Osh]. The general 
proof for Case (iii) (and hence Cases (i) and (ii)) was given recently by I.Cherednik 
[Ch]. 

If m = 1 then any Schrodinger operator is integrable by the definition. In two 
or more variables integrability is a very rare property. This is illustrated by the 
following result. 

Theorem 2.2. [OOS,OS] Let m > 2. Let L be an integrable Schrodinger operator 
defined by (2-1) with V = C Assume that the operators Li, 1 < i < m, are 
invariant under the symmetric group Sm+i acting irreducibly in C™", and their 
symbols generate the ring of Sm+isymmetric polynomials on C"^ (the operator Li 
is of order i-\-l). Then L coincides with (2-2) for the root system Am for some 
values of the parameters. 

Theorem 2.1 can be generalized to the matrix follows. 

Let be a complex simple Lie algebra, f| C be a Cartan subalgebra, C I)* 
be the root system of g. A""" be the set of positive roots. Ca, fa € fl be the root 
elements corresponding to the positive root a. Let U he a, representation of Q such 
that the space F = C/[0] of zero weight vectors in U is finite-dimensional. Define 
the matrix Schrodinger operator 

(2-3) Hg^u,u = A(, - ^ u{< X, a >)eafa, a; G 1), 

a6A+ 

where u is of type (i), (ii), or (iii) from Definition 2.3. Such operators are considered 
in [E],[EK1]. 

Theorem 2.3. The operator (2-3) is integrable. The symbols of its quantum inte- 
grals are generators of the algebra of Weyl group invariant polynomials on f). 

This result is proved in [E] for the special case g = slm+i, but the method used 
for the proof works for any Lie algebra. This method uses representation theory of 
the Lie algebra q in the trigonometric case, and representation theory of the affine 
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Lie algebra Q in the elliptic case. The quantum integrals of -ffg,i7,u are constructed 
from central elements of the universal enveloping algebra. We discuss this method 
in the Appendix. 

As a particular case, Theorem 2.3 includes Theorem 2.1 for the root system Am- 
Indeed, let us take g = slm and a special representation of g: 17^ — {zi...Zm+i)^C[^, 
fj, E C, with the action of g by linear transformations of variables (this represen- 
tation has no highest weight). All weight subspaces in are one-dimensional; 
in particular, V = U^[0] — C It is easy to compute that eafa\u[o] = A*(/" + !)• 
Therefore, if /x is chosen in such a way that = -|- 1), then operator (2-3) 
transforms into (2-2) for the root system A^- 

Krichever [Kr] introduced the notion of an algebraically integrable Schrodinger 
operator (see also [CVl]). Here we generalize this definition to the matrix case. 

Let L be an integrable matrix Schrodinger operator, let Li = L, ...,Lm be its 
quantum integrals, and let the symbols of Lj be piip/dxi., d/dxm)ld, where pi 
are algebraically independent polynomials. 

Definition 2.4. L is called algebraically integrable if there exists a matrix differ- 
ential operator Lq commuting with Li, ...,L^ with symbol po{d/dxi, ...,d/dxm)Id, 
such that for generic Ei, E^, £ C the polynomial po takes distinct values at the 
roots of the system of equations Pi{yi, y^) = Ei, 1 < i < m. 

For V = C this definition coincides with the one in [Kr],[CVl]. In the matrix 
case and m = 1 the property of algebraic integrability of differential operators was 
studied in [G]. 

It turns out that a Calogero-Sutherland operator is algebraically integrable for 
a discrete spectrum of values of the constants Cq. 

Theorem 2.4. [VSC] // Cq, = ^fi^ifJ-a + 1) < a, ct > for all roots a E A+, 
where //q is an integer depending only on the length of a, then the operator (2-2) 
is algebraically integrable for the rational and trigonometric potential. 

Conjecture 2.5. [CV2] Theorem 2.4 is true for the elliptic potential. 

Conjecture 2.5 is proved only for the case of the root system Ai. In this case, 
operator (2-2) is the Lame operator L = — Cp, and algebraic integrability of this 
operator is equivalent to the finite gap property, which takes place for C = //(// + 1), 
ji [Kr]. In this case, there is a quantum integral Lq of order 2// -|- 1. 

Now let us consider the case of the root system A^. Looking at the interpretation 
of the Calogero-Sutherland operator via the representation U^, we see why the 
integer values of jj, should be special: they are exactly those values for which the 
representation has a finite dimensional submodule or quotient module which is 
isomorphic to a symmetric power of C"^"'"-^ (or ((7"^+^)*). Since the zero weight 
vector is contained in this finite-dimensional module, we can use it instead of U^. 
Thus we observe that algebraic integrability occurs at those values of ^ where 
can be replaced by a highest weight module. This motivates the following general 
theorem which is the main result of this paper. 

Theorem 2.6. If U is a highest weight g-module then Hq^u,u is algebraically inte- 
grable for the rational and trigonometric potential. 

In Sections 3-6 we prove this theorem for the trigonometric case. The rational 
case can be obtained in the limit, so we don't discuss it. 



Note that Theorem 2.4 for the root system A^, is a special case of Theorem 2.6. 
Finally, we would like to formulate a natural conjecture concerning the elliptic 
case. 

Conjecture 2.7. Theorem 2.6 is true for the elliptic potential. 

This conjecture contains Conjecture 2.5 for the root system A^. We believe that 
it could be proved by applying the methods of this paper to the elliptic case and 
using the techniques of representation theory of affine Lie algebras and theory of 
vertex operators introduced in [E],[EK1]. 

3. Verma modules, Shapovalov form, intertwining operators. 

Let be a simple complex Lie algebra with triangular decomposition g = n_ © 
^ © n+. Fix an element A e f)*. Denote by M\ the Verma module over q with 
highest weight A, i.e. the module with one generator v\ and relations 

n+VA = 0, hv\ = (A, h)vx for h E i). 

We have the decomposition Mx = 0^<a-^a[a*] of M\ into the direct sum of 
finite dimensional weight subspaces Mx[iJ,]. Denote also = 0^<a-^a[a*]* the 
restricted dual module to Mx with the action of g defined by duality. is a 
lowest weight module with the lowest weight vector v'li^ of weight —A. 

We have a vector space decomposition U{g) = U{n-) ®U{l))®U {n+). Define the 
Harish-Chandra homomorphism : t/(0)[O] U{i)) by (j)\u{t)) — Id, and ^(fif) = 
if ^ e t/(£|)[0] can be represented as = giei for some gi e C/(£|),ei e n+. This in 
turn gives rise to a contravariant bilinear t/(f))-valued form T on ?7(n_) defined by 

^{91,92) = (f){uj{gi)g2), 

when gi,g2 belong to the same weight subspace of t/(n_), and J^{gi,g2) = oth- 
erwise. Here u is the Cartan antiautomorphism of Q defined by 

(^(ei) = /i, uj{fi) = Cj, uj{hi) = hi 

It is easy to see that this form is symmetric. 

As U{1)) can be identified with the space of all polynomials on [)*, we can intro- 
duce a symmetric contravariant C- valued form F on Mx defined by 

F{givx,g2Vx) = ^(9i,^2)(A). 

Let U be any weight g-module with finite dimensional zero weight space t/[0]. 
The completed tensor product Mx®U = Homc{A4^,U) has a natural g-module 
structure. We say that an element v ® u has order 77 if G Mx [A — ry] . Clearly, 
only elements of order 77 G may occur, where — '^+^ol- We say that 

v < rj if v ^ rj and rj — v & Q"*". 

Let u E U. Let • -^a Mx®U be an intertwining operator such that 
vx^ vx'^ u+{higher order terms} . 

It is clear from the intertwining property of $^ that u has to be a zero weight 
vector. 
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Proposition 3.1. 

(1) If Mx is irreducible then $^ exists and is unique for any u e C^[0]. 

(2) If Mx is reducible then $^ exists iff u & C/[0] is such that Ixu = 0, where 
Ix C U{n+) is the left ideal consisting of all g e U{n+) such that gv'^^ — 0. 

Proof. First consider the case when Mx is irreducible. Because Mx is freely gener- 
ated by vx over t/(n_), we only need to prove that the module Mx®U contains a 
unique singular vector of the form vx ® u+{higher order terms} . This is the same 
as to construct a map © : M^ — > U such that 0(f * = u and © is a n+-intertwiner. 
But M^ is a free ?7(n-|.)-module generated by vl^, so © can be uniquely extended 
from t;* ^ to the whole M^. Now let Mx be reducible. Then M^ is also reducible, 
and the universal enveloping algebra U{n^) contains the nonzero annihilating ideal 
Ix of the vector f 1^^. The same argument as above shows that Ixu = is a neces- 
sary and sufficient condition for the existense of the map ©. If it is satisfied we can 
put to be the singular vector of weight A in Mx<SiU, and then extend this 

map to the whole Mx- □ 

It is known that Mx is irreducible for generic A. For special A's Mx may be 
reducible, and it happens when the contravariant bilinear form on Mx is degenerate. 
Shapovalov [Sh] obtained an explicit formula for the determinant of this form: 

(3-1) det F^{X) = const Yl J] ((a, A + p) - ^(a, «))^('^-"«), 

aeA+ nGN 

where = {F/j). i, j = 1, 2, . . . , dimM^fA — //] is the matrix of the restriction 
of the form to Ma [A — ^u], K{ji) - the Kostant partition function, and the nonzero 
constant depends on the choice of basis in Mx[lJi]- 
Let 

(3-2) ^a(A) = (a,A + p)--(a,a). 

The conditions for reducibility of Mx can then be rewritten as 

Xn{>^) = for some a e A+,n e N 
Now we fix weight /x and let = max{n G N|i^(/x — na) ^ 0}. Denote 

(3-3) xm(a)= n n^n(A) 

aeA+ n=l 

We need the following 

Lemma 3.2. Matrix elements {F~^)^_. of the inverse matrix F~^ can be written 
in the form 

'''' Xm(A) 



for some suitable polynomials Pl^-{\). 

Proof. Shapovalov formula implies that matrix elements are rational functions in A 
with only possible poles in hyperplanes defined by X^(A) = 0. Our goal is to show 
that only simple poles may occur. 
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Fix a e A+, n < tt-q,. Take A such that Xm('^) = iff m = n, /? = a. Then M\ is 
reducible and contains a unique maximal submodule of M\ , generated by a singular 

vector V\-na- 

Fix 2; G [)* such that (a, 2;) ^ for any a G A""", and let t be an independent 
variable. Using the [/([)) -valued bilinear form T we can introduce a new C[t]-valued 
bilinear form F* on M\ defined by 

F\givx,g2vx) = ^{gi, g2){^ + tz), gi,g2 G f/(n_) 

Clearly, specialization 1 1— > gives the usual Shapovalov form. 

Denote N = dim Mx [A - //] = K{p,) , M = dim [A - ^f] = (/U - na) . Choose a 
basis Vk in Mx[X — /j] so that {t'i}, z = 1, 2, . . . , M, would form a basis for M^[A — //]. 
Then the matrix elements {F^) . . will be divisible by t if i < M or j < M : 





t ■ fl,M 


t ■ /l,M-|-l 


■ t-fl,N\ 


t ■ /m,i 


t ■ fM,M 


t ■ fM,M+l ■ 


t ■ fM,N 


t ■ /m+1,1 • 


t ■ fM+l,M 


fM+l,M+l ■ 


fM+l,N 


V t-fN,l . 


t ■ fN,M 


fN,M+l 


fN,N / 



where fij are some polynomials in A, t. It is clear now that the determinant of 
any (A^ — 1) x (A^ — 1) submatrix of Fj^ is divisible by t^~^. Shapovalov formula 
implies that det is divisible by exactly Mth power of t, which means that when 
we compute the matrix elements of (F^) ^, only simple poles will be allowed when 
t = 0, or, equivalently, {F~^). ., will have at most simple poles on the hyperplanes 
X"(A) = 0. Repeating this argument for all n, a we prove the lemma. □ 

We can apply this result to get more information about the intertwining operator 
In the proof of Proposition we defined $^va as a map : — > U. We would 
like to obtain a more explicit formula for ^^vx as an element of Mx®U. 

For any basis Vk = g^'^x, g^ ^ U{x\.-) of Ma[A— //] we have the basis of M^[— A-|-//] 
given by f ^ = (a;5f^)f * where uj is the Cartan involution. It is clear that (f *, Vk) = 
F{vi,Vk).Oz 

Introduce another basis Wk which is dual to in the usual sense, i.e. (f *, Wj) = 
Sij. These two bases Vk and Wk are related via the matrix: 

i 

i 

It is clear that in this notation 

= vx®u+- ■ ■+J2wk®i^g'^)u+- ■ ■ = vx®u+- ■ {F-^)ki9k^^®i^9i)u+. • • 
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Corollary 3.3. Suppose we have an intertwiner : Mx — > M\®U, where U is a 
highest weight Q-module with the highest weight 9. 

(1) There are no order ^ terms in the expression for unless <9. 

(2) If n < 9 then the order jj, part of a can be written as 

(3-4) Y.^F-^)ki-9k^^®^^9t)u. 

k,l 

where 

^ ^ Xm(A) 

for some polynomials P^i{X). 

(3) The condition of existence of can be written as Xoi^) 0- 

(4) for any basis g'j^ o/C/(n_) we can choose polynomials Sj^i{X) so that 

(3-5) nvx = E E ?7iY^^^ ® 



For a rational function R, represented as a ratio of two polynomials R = ^ we 
set deg R — deg P — deg Q. 

Note that all coefficients (F~^)^^ are of negative degree in A. Later we will work 
with A in the hyperplanes {a, A) = const, so we introduce notation 

(3-6) Aq. = \ ' l a, A"*- = A - Aq,, 

[a, a) 

so that A^ is a (dim[) — 1) -dimensional vector and {a. A"*-) = 0. 
We will use the following 

Proposition 3.4. When restricted to the hyperplane {a, A) = C, matrix elements 
(F^^)^^ are rational functions in A^ of nonpositive degree, and only constants may 
occur as terms of degree 0. 

Proof. 

We choose a special basis in C/(n_ )[//]. For any sequence uj of positive roots 

/3i > > • ■ ■ > /9r5 where > denotes now the lexicographical order, such that 
=11, set = fp^ . . . fp^. Set degX^ = Card{{k\Pk «})• We also write 
fp = g-p, ep = gp for a positive root 

The set of X^^'s is a basis in C/(n_)[/x]. We also have 



(3-7) deg^(X^,,X,,)(A)< 



degX^i +degX^^ 



Indeed, we can only raise the degree by commuting some ep and fp for P ^ a, 
which results in the term (/?, A) -|- const, which is linear in A"*-. Note also that 
commuting with e^ or f^ will not increase the total number of terms ep and fp for 
all (3 ^ a. Therefore, the maximal degree cannot be greater than half the original 
number of terms ep and fp, (3 ^ a. This proves formula (3-7). 



The determinant of the form in the hyperplane (A, a) = const is equal to 
det F^(A) = const J] H ((/^' A + p) - 

where the constant depends on C. Then 

iV= J]^K(^-n/3) 

is the degree of the determinant as polynomial in A"*-. By the same argument as in 
[Sh], from (3-7) it follows that the A-'"-dcgree of any minor of the Shapovalov matrix 
cannot exceed A^. Moreover, commuting with Ca or does not change the set of (3 
mod a, and therefore any term of degree exactly N has highest term proportional 
to that of the determinant, which proves the Proposition. □ 

4. Matrix Trace, ^-function and its properties. 

Fix a highest weight g-module U with highest weight 6 and finite dimensional 
zero weight space U[0]. Consider a new operator 

(4-1) ^^ = X^(A)^^ 

From Corollary 3.3 it follows that 

M k,l 

This expression allows us to define $^ even for A such that Xe(A) = where 
itself was not defined. It is clear that $^ is an intertwining operator for all A, and 
it has the property 

(4-2) ^a'^a = vx ® X6((A)w -|- {higher order terms} . 

Let ^{X,x) be a £'n(i(C/[0])-valued function on ^* x ^ defined by 

(4-3) =^5-1*1!) 

Proposition 4.1. The ^-function defined above has the following properties: 
(1) 

(4-4) ^'(A,x) = e<^+^'">P(A,x), 

where P{X,x) is a End{U[0]) -valued polynomial in X with the highest term 

n («,A)"«.M 

and we put for brevity na = n^. 
(2) If {a, X-\- p) — f (a, a) = for some a e A'^ ,n = 1, 2, . . . , tt-q, then 

(4-5) '^{X,x) = ^{X-na)B^{X), 

for some -B^(A) e End{U[0]), which is rational in X of nonpositive degree, 
and only constant operators may appear in it as degree zero terms. 
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Proof. The first part is clear from the formula 

and the fact that all the Tr\Mxl\-i3]i^\) some combinations of 5'^;'s, and there- 
fore polynomials in A. Their highest terms are obviously all equal to naeA+ ('^' ^)"^°' ■> 
so the highest term of P(A, x) is equal to 

We now prove the second property of the \l/-function. Let (a, A + p) — ^(o;, cc) = 
for some a e A+, 1 < n < ria, but A + p) — /5) 7^ unless P = a,m = n. 
Set = ncu. 

From Corollary 3.3 it follows that a has no order u terms unless u > na. In 
particular, there are no order zero terms. On the other hand, ^^vx has to be a 
singular vector. Therefore we must have 

(4-6) — v\—na 'S> u -\- {higher order terms} , 

where v\-na is the unique singular vector generating the submodule — M\-na- 

This implies that ^\ is a triangular operator: ^\M\ C M\ ® U, so 

Tr\Mx{ne'') = Tr\Mi{ne''). 

Let ^'xVx-noi = v\-na ® w + {higher order terms} . Using the fact that = 
Mx-na we see that 

(4-7) TrlMxinen = Tr|M,_„„ (^^-^ae^) = ' ^ " 

Xe{X - na) 

It is an easy calculation to show that 

Introduce a linear operator 

(4-8) B^iX) = Y: ■ 9M9n e EndiUm 

^ Xe[A - na) 

We can rewrite (4-7) as 

^'(A,x) = ^{X-na,x)B^{X), 

To complete the proof we only need to show that -B^(A) satisfies the required 
condition. It is clear that -B^ ( A) is rational in A and is not singular in the hyperplane 
{a, X + p) — ^{a, a) = 0. As we can rewrite 

the rest follows from the Proposition 3.4. □ 

Now we can introduce our main object of study. Set k = na- Put 

(4-9) iPiX,x) = 2-^i^-p,2x) 

The properties of ^^-function can now be rewritten in the following form: 
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Corollary 4.2. 

(1) ijj-function can be represented as 

(4-10) tp{X,x) = e<^'''^P{X,x), 

where P{X, x) is a polynomial in X of the form 

P{X,x) = Y\ (cK, A)"" + {lower degree terms} . 

(2) // {a, A) = for some a e A"*", then for s = 1,2, ... ,na. 

(4-11) i/;{X + sa,x) = ipiX- sa,x) ■ B^{X), 

where B^{X) = B^{^-p) is a rational End{U[0]) -valued function of X and 
can be represented as 

B"{X) = bg + {lower degree terms} 

for some constant bg e End{U[0]). 

5. Uniqueness of the V'-function. 

In this section we prove the uniqueness property of the function '0(A,x), satis- 
fying (4-10) and (4-11). 

Proposition 5.1. 

Suppose we have an End{U[0]) -valued function 

0(A,a;) = e<^'^>g(A,a;), 

where Q{X,x) is a polynomial in X, satisfying (4-11)- Then the highest term of 
Q{X,x) is divisible by 

n («,A)"". 

Proof. Consider the highest term of Q{X,x). We need to show that it is divisible 
by {a. A)"" for any a e A+. 

Fix an q; e A""". We can uniquely represent Q{X,x) as 

L Ki 

(5-1) QiX,x) = Y,Y.^-Qki{X^,^), 

where QkiiX-^,x) are homogeneous £'n(i(C/[0])- valued polynomials in A"*- of degree 
I. 

The highest term of Q{X, x) will be some combination of the terms X^Qkii{X, x). 
We claim that it is enough to show that Kl > ria- Indeed, it will follow then that 
the highest term will have degree at least L-\-na, and therefore all terms of the form 
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X^Qkii{^,x) contributing to the highest term must have Ki > L + Ua — I > n, 
which proves the statement. 

By our assumption 'iIj{X, x) satisfies (4-11), so we can write 



/ L Ki 
\l=0 k=0 

(5-2) e-^<"'^> (j2 Yl'^-s)''Qkl{^^,x) \ (6^ + {lower degree terms} ), 

\l=0 k=0 / 

where {lower degree terms} are understood with respect to A"*-. 

We can consider homogeneous parts of (5-2) of degree L in A"*-. FormaUy, given 
a function /(A-*-), we consider 

lim p — . 

t^oo t^ 

This gives us 




(5-3) e^<-'-> (f]s'^gfcL(A^,a;) ) =e-^<"'=^> ( E(-s)'^QfcL(A^,a;) ) 6« 

for s = 1, . . . , Kl. 

The rest is based on the following 

Lemma 5.2. Consider a homogeneous system of N linear equations on K vector 
variables Ak{z) e , which are meromorphic in some additional parameter z : 

K 

(5-4) J]s^(e- + C,)Afe = 0, 

k=l 

s = l,...,N,Cs e MatM{C). 

If K < N then this system has only trivial solution Ak{z) — 0. 

Proof of lemma. 

We can think of this system as a system of linear equations on KM variables 
{Ak)m and rewrite (5-4) in the block- matrix form 



/ (e^-Id + Ci) {e^-Id + Ci) 

2(e2^-/d+C2) 22(e2^ -Id+Ca) 

\N{e^' ■ M+Cn) N^{e^' -M+Cn) 



(e^-Id + Ci) \ /Ai\ 



2^(e2« -Id+Cz) 



N^ie^' -M+Cn)/ \Ak) 



Suppose K < N. Then the determinant of the submatrix, consisting of first K 
blocks (or, equivalently, first KM equations) is an entire MatM (C)-valued function 
of z with the asymptotics as 2; — ^ -|-oo 
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det 



/ (e^ ■ Id) 
2{e^' ■ Id) 



(e^ ■ Id) 
2^e^'-Id) 



\K{e^'-Id) K^{e^^-Id) 



(e^ ■ Id) \ 
2^{e^'-Id) 

K^{e^'-Id)/ 



which is equal to 



where the nonzero constant 



r K(K+1) . 

D ■ e 2 



M 



D = det 



/I 

2 



1 

22 



is the Vandermonde determinant. 

Therefore, this determinant is a nonzero entire function, which imphes that for 
generic z it is not zero, so the system has only trivial solution. The meromorphic 
functions Ak{z) are equal to zero for generic z, and therefore must be identically 
equal to zero. 

The lemma is proved. □ 

We now apply Lemma 5.2 to the system given by (5-3), ior N = Ua, K = K^ + l, 
M = dimC7[0] and setting Cs = (-1)''+^ {b"^)* , z = = ^^Oi, where is the 
transposed matrix A. 

Consider the rows of the matrix, corresponding to '0(A-'-, x), and transpose them 
so that they become columns. By (5-3) all these columns satisfy the system of 
equations (5-4), and as ■^-function is not identically equal to zero, it implies that 
the system (5-4) has a nontrivial solution. By Lemma 5.2, we have Kl + 1 > Ua, 
or, equivalently, Kl > rio,. 

The proposition is proved. □ 

Corollary 5.3. 

Any End{U[0]) -valued function 

</.(A,a;) = e<^'^)Q(A,a;), 
satisfying (4.-11), can be represented as 

(j){\x) = q{X)ip{X,x) 
for some End{U[0])-valued polynomial q{X)- 

Proof. The proof is similar to the proof of the Lemma in Section 1 of [CVl]. We 
use induction on the degree of Q{X,x). 
If deg(5(A,a;) < degP(A,a;), where 



V'(A,a;) = e<^'^>P(A,a;), 



then by proposition we have Q(A, x) = 0, so we can take g(A) = 0. 
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Suppose we have proved the statement for aU polynomials of degree less than 
the degree of Q{X,x). By Proposition we can find a £'n(i(C/[0])- valued polynomial 
gi(A) such that 

highest term of Q{\) = gi(A) ]^(q;, A)"'" 

a 

Consider the function 

01 (A, = 0(A,a;) - qi{\x)ii{\x). 

Obviously, it satisfies (4-11). Moreover, it can be represented as 

(/>i(A,a;) = e<^'">g(A,x), 

where polynomial Q{X, x) has degree smaller than that of Q{\, x). 

By induction hypothesis we can introduce a £^n(i(L''[0])-valued polynomial 52 (A) 
such that 

0i(A,x) = q2{\)il){\x). 

The polynomial q{\) ~ gi(A) + g2(A) satisfies the required property. □ 

Corollary 5.4. The function iIj{X, x) , satisfying both (4-10) and (4.-11), exists and 
is unique. 

Proof. It is a direct consequence of (4-10) and Corollary 5.3. □ 

6. Existence of differential operators 

The properties of the -i/'-function obtained in Chapters 4,5 are very close to the 
axioms in [CVl]. The function satisfying these axioms was used to construct a ring 
of differential operators that contained dim f) algebraically independent operators, 
corresponding to the generators of the ring of W-invariant polynomials, but was 
bigger than the ring generated by those operators. 

Here we apply these ideas to construct a similar ring of matrix differential oper- 
ators and thus prove Theorem 2.6. 

Theorem 6.1. For any End{U[0]) -valued polynomial Q{X) satisfying the property 



(6-1) Q{X + so) = Q{X — sa), s = l,2,...,na 

whenever (a. A) = 0, there exists a differential operator Dq with coefficients in 
End{U[0]), such that 

(6-2) DqiP{X,x) = Q{X)iP{X,x). 

The correspondence Q{X) Dq is a homomorphism of rings. In particular, all 
Dq commute with each other. 

Proof. 

We use induction on the degree of Q{X). If degQ(A) = 0, then Q{X) — const, so 
the operator Dq will be just the operator of multiplication by this constant. 
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Suppose we have proved the theorem for aU polynomials of degree less than that 
of Q{X). Let the highest term of Q{X) be equal to 

highest term of ((5(A)) = ^a(n)A("), 

(n) 

where (n) is a multiindex, a(n) e End{U[0]). Consider the operator Dq defined by 

~ V- 

(n) 

It has the property that 

Dqi(;{X,x) = Q{X)i(^{X,x) + {lower degree terms} , 

and it also satisfies (6-1). Consider the difference 

(/.(A, x) = DqiP{X, x) - Q(A)V'(A, x). 

It satisfies (4-11) and therefore can be represented as 

(l){X,x)^Q{X)ij{X,x) 

for some i??i(i([/[0])-valued polynomial Q{X) such that degQ{X) < degQ{X). By 
induction hypothesis we can introduce an operator Dq such that 

(j){X,x) ^ DqiP{X,x) 

and the operator 

Dq = Dq + Dq 

has the required property, which completes the proof of the induction step. 

The assertion that the constructed correspondence is a homomorphism of rings 
follows from the fact that the operator Dq^q^ —Dq^Dq^ annihilates the '(/'-function 
for any A, and therefore has to be identically zero. □ 

Among the polynomials (5(A), satisfying (6-1), are all M^-invariant polynomi- 
als pi(A),...,Pr(A). It is known that they are algebraically independent, and the 
ring generated by corresponding differential operators is a ring of polynomials in 
generators Dp,, . . . , Dp^. 

There are also other polynomials, satisfying (6-1), which are not VF-invariant. 
They give rise to differential operators which are not VF-invariant and therefore do 
not belong to the ring generated by Dp^ , • • • , Dp^. 

In particular, all polynomials contained in the ideal generated by 

Qo(A)= n n((«'^)'-^' 

q:€A+ n—l 

satisfy the (6-1). 
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Proposition 6.2. The differential operator corresponding to the invariant polyno- 
mial pi{X) = (A, A), is equal to 



This fact can be proved by a direct computation. Another proof of it using the 
relationship between the center of U{g) and commuting differential operators is 
sketched in the Appendix (see also [E]). 

The operator (6-3) coincides with the generalized Calogero-Sutherland operator 
(2-3) for the trigonometric case. We have shown that this operator is algebraically 
integrable. Theorem 2.6 is proved. 

Appendix 

In conclusion we briefly describe how to construct quantum integrals of Hg^jj^u 
with trigonometric potential from central elements of U{q). This construction works 
for arbitrary module U, not necessarily highest weight. 

Proposition Al. [E] Let X G U{g) be an element of degree 0, i.e. [h,X] = 0, /i e 
[). Then there exists a unique matrix differential operator 'D{X) with End{U[0])- 
valued coefficients such that for any A e f)* and any intertwining operator $ : Mx — > 
Mx^U 

Tr\M,{^Xe^) = V{X)Tr\Md^e^). 

The proof of this theorem and a recursive construction oiT>{X) is given in [E]. 
We illustrate the idea of this construction by computing D{EF) for g = SI2 {E, F, H 
are standard generators of g). The main trick is to carry E around the trace, using 
the intertwining property of $ and the cyclic property of the trace: 



Tr|M,($£^i^e") = Tr|M,((£^ ® l)$Fe") + ETr\M,{<^Fe'') = 
Tr\M^{'^Fe''E) + ETrlM^'^Fe'') = e<'^^''>Tr\M^{^FEe'') + ETr\M^{,^Fe'') = 
e<^'^>Tr\M^{^{EF + F)e^) + ETrlM^i^Fe"") = 

(Al) e<"'^>Tr|M;,(^£^i^e^) + e<"'^>^Tr|M^(*e^) + STr|M;,(^i^e^), 
(where a is the positive root of g) which implies 

(A2) Tr\MA^EFe-) = ^ _ (^e<"-">^rr|M,(^e") + ^Tr|M,(^Fe")^. 

Further, we have 
(A3) 

Tr\M,{^Fe^) = Tr\M,{{Fm)<^e^)+FTr\Md^e^) = e-<"'">Tr|M,($Fe^)+FTr|M,($e^), 



which implies 

(A4) TrlMd^Fe-) = ^ _ ^-<^,.> ^^rlM, (^e^)- 
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Combining (A2) and (A4), we find 



Thus 

g<a,a;> ^ 

(A6) - YZr^^^Q^ + (l - e<«'^>)(l - e-<«'^>)^^' 



In general, it is not true that V{XiX2) equals either V{Xi)V{X2) or P(X2)P(Xi). 
However: 

Proposition A2. [E] IfXi belongs to the center ofU{g), then for any X2 one has 

V{X2)V{Xi)=V{XiX2). 

This proposition follows from the fact that if $ is an intertwining operator then 
$Xi is also an intertwining operator. 

Corollary. If Xi,X2 are both in the center ofU^Q), then V[Xi) andV{X2) com- 
mute with each other. 

Let D{X) be the differential operator obtained from "D^X) by conjugation by the 
function Tr|Af_,(e-), i.e. defined by D{X)^{x) = Tr|M.,(e^)P(X)(^^^^g^). 
The following statement is checked by a direct computation: 

Proposition A3. [E] Let B be an orthonormal basis of q, and Ci — "^aeB*^^ 
the Casimir element. Then 

(A7) D(Ci) = A. - y T^T^^ ^ const. 

aeA+ ^ ' ' ' 



Proposition A4. Let Ci,...,Cm be algebraically independent generators of the 
center of U{q). Let Li = D{Ci). Then Li,...,Lm are algebraically independent 
quantum integrals of the Schr — "odinger operator Li given by (A7). The symbols 
of Li generate the algebra of Weyl group invariant polynomials on [). The function 
\1/(A, x) is a joint eigenfunction for Li, L^, with eigenvalues (f){Ci){X), where cj) 
is the Harish- Chandra homomorphism defined in Section 3. 

Observe that operator (A7) transforms into jHg^u,u + const when one makes a 
change of variables x ^ 2x. Therefore, we have 

Proposition A5. Let Li be the operators obtained from Li by replacing x with 2x. 
Then Li — ^Hg^u^u + const, and Li, ...,Lm are algebraically independent quantum 
integrals ofHg^u,u for trigonometric u. The function '0(A, x) is a joint eigenfunction 
of Li, Ljn. 

This implies Theorem 2.3 in the trigonometric case. 

Finally, we observe that if C is a central element of U{g) then by Propositions 
A4,A5 one has D{C) = Dp, where Dp is defined in Section 6, and p is the Weyl 
group invariant polynomial given by p{X) — 0(C)(2(A + p). This proves Proposition 
6.2. 
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